Understanding spin Hall effects from the motion in SU(2)xU(l) fields 
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We derive the classical counterpart of a previously obtained quantum mechanical covariant "con- 
tinuitylike" equation for the spin density, and present an intuitive picture for elucidating the non- 
conservation of the spin current. This reveals the equations of motion for a particle with spin under 
the Yang-Mills field (in certain semiconductors) as well as the Maxwell field, from which the condi- 
tion for an infinite spin relaxation time is drawn out directly. As a concrete example, we discuss the 
procession of the spin orientation in spin Hall effect with the so called ReD field, which undergoes 
a circle with the frequency dependent on both the strength of the spin-orbit coupling and the ini- 
tial velocity. The anti-commutation of the Pauli matrices is found to be crucial in simplifying the 
equations of motion in the view of quantum mechanism of the same topics. 

PACS numbers: 72.25.Dc, 72.25.-b, 03.65.-w, 85.75.-d 



The spin Hall effect 0, Q is regarded to be practically 
useful for the rapidly developing field of spintronics , in 
which the spin manipulation, polarization as well as the 
detection play crucial roles. Theoretically, the spin Hall 
effect was predicted for p-type semiconductors |4j where 
the up and down spin particles drift in opposite directions 
due to an "effective magnetic field" originating from the 
Berry phase curvature, and a constant spin Hall conduc- 
tivity was speculated for n-type semiconductors |5| with 
the help of either the Bloch equations or the Kubo for- 
mula. Although the absence of the Hall voltage makes 
it difficult to detect the pure spin Hall effect directly, 
the spin accumulation in nonmagnetic semiconductors 
has been observed experimentally by two groups @, 
recently. 

The nonconservation 0, of the spin density in the 
presence of spin-orbit couplings brings about some diffi- 
culties relevant to the theoretical analyses |^ llfllri| . such 
as the definition of the spin current. A covariant form for 
the continuitylike equation for spin density and current 
was given 9] in the terminology of Yang-Mills gauge po- 
tentials, and it was shown [T^ to paly an essential role in 
guaranteeing the consistency of a generalized Kubo for- 
mula for the linear response to non-Abelian fields. Since 
the work mentioned above are all based on quantum me- 
chanics, one may ask what is the classical counterpart of 
the covariant form of the continuitylike equation for the 
spin density, and what is the classical interpretation for 
the non-conservation of spin current? 

For a semiclassical understanding of the spin Hall ef- 
fects in a two-dimensional electron system with Rashba 
coupling, an external electric field was shown to exert a 
transverse force 0, IT3L ITU on a moving spin, where the 
calculation in Ref. 14] is done on the basis of Ehren- 
fest principle and the discussion in Ref. |l3l | is carried 
out with the help of an intuitive picture that the center 
of a wave packet gets deflection in the transverse direc- 
tion. The general form of the forces acting on spin and 
spin current for a variety of systems considered in cur- 
rent literatures was given in Ref. [t| . The time-dependent 



Rashba spin-orbit coupling was ever proposed to create a 
force acting on opposite spins in opposite directions 0] . 
This makes it possible to produce the Yang-Mills field in 
certain semiconductors according to the viewpoint of our 
formulism [9j . It is natural to ask whether our formulism 
has advantage in making sense of the roles that spin and 
charge play in spin Hall effect. It is therefore an emergent 
task to elucidate the aforementioned issues via a unified 
picture. 

In this letter, we investigate the motion of a particle 
with spin in the Yang-Mills field from a classical point of 
view. We start with a revisit to the classical counterpart 
of the continuitylike equation for the spin-current density. 
Such a topic was ever investigated by various authors |l6( 
but has not been correctly exposed yet. We derive the 
same form as ever proposed quantum-mechanically in the 
previous paper [i| , which presents a clear picture of the 
non-conservation of the spin current. On the basis of this 
result, we obtain the equations of motion for a particle 
with spin in the presence of the Yang-Mills field as well 
as the Maxwell field. A straightforward condition for 
an infinite spin relaxation time is drawn out. Then we 
discuss the spin Hall effect as a concrete example and 
obtain the trajectory of the spin for the equal strength 
of Rashba and Dresselhaus couplings. We also discuss 
the quantum counterpart of the equations of motion and 
indicate that the anti-commutation relation of the Pauli- 
matrices makes the problem easy to solve. 

We start with considering a moving top (classical anal- 
ogy of spin) which rotates at a certain rate. A continuum 
constituted by such kind of tops is completely charac- 
terized by a local velocity field v(r, t), a local particle- 
density field p(r, t) together with a local alignment field 
N(r,t) [13. As illustrated in Fig. QJ, the time evolu- 
tion of N(r, t) is determined by comparing N at different 
times at the same place, i.e., N(r,t + At) — N(r,t) = 
Q x N(r,t)At, while its spatial deviation is determined 
by comparing TV at different places simultaneously, i.e., 
N(r + Axi,t)-N(r,t) = f2j x Y(r,t)Ax l . Hereafter, the 
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FIG. 1: (color on line) The left scheme depicts the time- 
evolution of iV at a point r; the right scheme illustrates the 
spatial deviation of N's by comparing the fields at two neigh- 
borhood points r and r+Axi for which the parallel translation 
is inevitable. 



Latin indices run from 1 to 3, and the repeated indices 
are summed over. The vector fields uj and are natural 
consequences of N being a unit vector, then we have 
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— JV(r,t)=£J(r,t) x N(r,t), 
^-iV(r,t) = O i (r,t) x N(v,t). 



(1) 



By making use of these two relations together with the 
, .dp dj, 

density conservation — + — — = 0, wc obtain a Continu- 
ed OXi 

itylike equation 
, 9 



)s + (- — n t x)ji = o, 

ox i 



as long as the nature definitions of spin density a — p N 
and spin-current density Ji — pvi N are employed. Com- 
paring with the quantum mechanical results [9( , one can 
recognize that cu and Qj correspond to the Yang-Mills 
gauge potentials r/Ao and — rjAi, respectively, which have 
been shown |fj| to describe the spin-orbit coupling such 
as Rashba Dresselhaus 0] coupling and etc. As 
we are aware, a clear and direct physical meaning of the 
non-conservation of spin current has not been exposed 
before. 

The above results are suggestive for us to elucidate the 
equations of motion for a charged vector ft of constant 
magnitude (classical analogy of a charged particle with 
spin) moving in the space undergoing both the Maxwell 
and Yang-Mill fields 



7]{Ao - v.A) x n(t), 



dn(t) 
dt 

dv ■ — * -* 

= & ■ + eE i + e ijk v j(Bk ■ n{t) + eB k ) 



(3) 



The first equation for the spin orientation is obtained 
from Eq.JU by adopting N(r,t) — n(t)5(r — r(i)) with 
f (t) being the trajectory of a single particle. The sec- 
ond equation is due to the fact that the translational 
motion is governed by both the Lorentz force caused by 
the Maxwell fields Ei and Bi and the force by the Yang- 
Mills fields £ i and Bi which are matrix- valued vectors [3 . 



Hereafter, n(t) is written as n for simplicity. To avoid 
ambiguity, we need to set up a picture with two distinct 
spaces. One is the SU(2) Lie algebra space (we call it spin 
space hereafter) in which the Yang-Mills gauge poten- 
tials and fields are defined, the other is the conventional 
spatial space. The "coordinate bases" of the former are 
{f 1 ,f 2 ,f 3 } with 2r referring to the Pauli matrices, e.g., 
Bj = Bjr 1 + Bfr 2 + Bfr 3 : while those of the later are 
{e x ,e y ,e z }i e.g., v = v x e x + v y e y + v z e z . In the lan- 
guage of the gauge potentials, the Yang-Mills "electric" 
and "magnetic" fields can be expressed as 
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-d Ai - diAo + r)Ao x A l 



Bi = e i: jkdjAk + ^ijkAj x A k - 



(4) 



It is worthwhile to note that the non-Abelian fields can 
be non- vanishing even when the gauge potentials are con- 
stant. 

In the view of the analytical mechanics, the equations 
of motion @ can be formulated in the Hamiltonian for- 
mulism with H = (Pi ~ eAi — rjAi ■ n) 2 + cAq + j]Aq ■ n. 
As a standard procedure, one needs to choose two sets of 
canonical coordinates and their conjugations (canonical 
momentums), namely, (cp,n z ) with cf> = tan _1 (— n x /n y ) 
and (ri,Vi). One pair of the canonical equations can be 
easily derived 



d(j) 
~dt 
dn z 
~dt 



dH 



® ? = 7T~ - ■\(( A o- v i A i) n2 -n z (Ao-ViAi)-n), 



dn 



dH 



with n 2 , = 



(r)(Ao - ViAi) x n) , 



n 2 , which is exactly the first equation of 



Eq. © after some algebra. The other pair of the canon- 
ical equation respect to (r^,^) can be shown to be the 
second one of Eq. ©. 

Some immediate consequences can be obtained from 
Eq. ©. (i) The first equation for the time rate of n 
clearly manifests that the vector n does not process when 
it is parallel to Aq— ViAi, or more specially, Aq— ViAi — 0. 
In such cases, the spin orientation keeps unchanged, 
which results in an infinite spin relaxation time. A typical 
example is that the infinite spin relaxation time occurs 
in the ±[1, ±1, 0] direction when Aq — and a = ±/J (a 
and (3 refer to Rashba and Dresselhaus coupling strength, 
respectively, in the quantum case), one of which was dis- 
cussed in Ref. [2(|. This is analogous to the case in the 
classical electrodynamics where an electron moving in the 
uniform orthogonal electromagnetic fields with velocity 
v = E x H/B 2 does not feel the Lorentz force, (ii) In 
virtue of the coupling between n and the Yang-Mills field 
in the second equation of © , the time rate of n leads to 
the time-dependent effective fields even when the Yang- 
Mills fields are time- independent, which is quite different 
from the motion in the Maxwell field. 

For many real physical systems, the gauge potentials 
Ai are constant. If we focus on the contribution of the 
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Yang-Mills field and set the Maxwell one to be zero, the 
second equation of © reduces to 



i— — 

dt 



-rjAi ■ (Ao — v jAj) x n, 



which gives rise to a relation between n and Vi 

1 



Vi(t) 



-At ■ n(i) + C h 



(5) 



(6) 



where Cj = VQi + -^Ai ■ uq are determined by the initial 
values i>j(0) = t>oi and n(0) = {tiq x , no v , uq z ). Conse- 
quently, one only needs to solve one equation, 

dff —* —* Tt -* —* 

— = r)(Ao - CiAi) x n + —(Ai x n)(A l ■ n). (7) 

dt m 

As a concrete example, we consider a two-dimensional 
system with uniform Yang-Mill gauge potentials 9] 



r 2m, ? 2m 
A = — /3, a, , .4^ = r 



(a, (3,0), A o = 0, (8) 



which correspond to Rashba and Dresselhaus couplings 
in certain semiconductors. Accordingly, the equations of 
motion are explicitly written as 



dn x 
~dT 



2m 

? 

V 



driy 2m 
dt T] 
dn 7 2m 



dt 



n 



2m 

) 

V 



ad - f3C 2 
/3d - aC 2 - 
ad - pd - 
Ci/3 — C 2 a 



4a/3n x + 2(a 2 +/3 2 )n y 



Aa(3n y + 2{a 2 +[3 2 )n x 



4a(3n x 



Tf 

4a/3n y 



(9) 



For a = j3 which is called ReD field in Ref. |2£|, we can 
solve these equations analytically: 

n x = -asin(wi + <p) + -(n 0x + n 0y ), 
n y = asm(u)t + (p) + -(n x+noy), 



= v / 2acos(w< + ip), 



(10) 



where a = \J n^ z /2 + (tiq x — riQ y ) 2 /4 and tan tp = (no y — 
no x ) / (\/2 noz) are determined by the initial conditions. 
It is clear that the tip of n experiences a cyclotron ro- 
tation with frequency u — 2y/ | ma (v ax ~ v 0y ), which is 
shown in Fig. J5J for the initial condition Hq = (1,0,0). 
The instantaneous velocity solved from Eq. JJjJ is just its 
initial value v x = vq x and v y = vq v , i.e., the electron 
undergoes a motion with uniform velocity. This is due to 
the time-dependent parts of n x and n y only differ from 
each other by a minus sign. Specially, when v x = v y , 
the spin vector n does not process since u> = 0, which 
recovers the result in Ref 
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FIG. 2: (color on line) The rotation of n is depicted with the 
initial conditions no = (1,0,0), where the x-y projection is 
also shown to elucidate this motion. 



When an external electric field E — (E x ,E y ) is ap- 
plied, which mimics to the usual spin Hall effect in cur- 
rent literature, we obtain Vi(t) = —^Ai-ri(t) + Ci + -^Eit. 
The equation of motion for n is almost the same as Eq. (JJJ 
if substituting d by d(t) = d + —Eft, For the electric 
field is sufficiently weak, the perturbation theory is ap- 
plicable and we can expand the spin unit vector n in the 
power of the electric field, 



ft = n<°) 



(11) 



The Oth-order results take the form of Eq. I|10f) . while the 
linear corrections are 



(i) _ 



n 



-\{E) t 2 cos(ujt + ip), 
X(E) t 2 cos(uit + ip), 
-y/2X(E) t 2 sin(ujt + ip), 



(12) 



where X(E) = aea{E x — E y )^/2/r]. The form of the entire 
motion for ft is much similar to the case without exter- 
nal electric fields, except for a time-dependent amplitude 
a(t) = [a 2 +A 2 i 4 ] 1/2 and a phase (p(t) = ip+taxr 1 {Xt 2 /a). 
Even though the rotation of ft is a bit of intricacy, its 
translational motion only involves one more term linearly 
dependent on time, 



''Ox 



V y 



— E x t, 
m 



(13) 



Furthermore, we turn to the quantum picture of the 
motion of a single particle with spin in the background of 
the Yang-Mills field with the Hamiltonian H = ^(pi — 

eAi — 7]A\ ■ f) 2 + cAq + r/Ao ■ f 9]. The time rate of 

its dynamical momentum 7^ = pi — eAi — rjAi ■ t is 
given by the Heisenberg equation of motion dfci j dt = 
i [%i , H] = eEi + Tj£i-f+ \e ijk {vj , eB k + r\B k ■ f }, where 
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the curve parentheses denote the anti-commutator. Sim- 
ilarly, the equation of motion for the spin operator reads 
df/dt = jf-[r,H] = r)Ao x r— %{vj,Aj x r}. Obviously, 
these equations are the quantum counterpart of Eq. Q 
with a r\T to n correspondence. In the previous classical 
case, it is difficult to solve the equations of motion since 
they are coupled to each other. However, things become 
easier in the quantum case. Substituting Vj = Ttj/m — 

— (pj — eAj — rjAj ■ t) into the equation of motion for r, 
we obtain df / dt — r/Ao x r — %{pj — eAj , Aj x r} in which 
the anti-commutation relation {r a ,T fc } = ^S ab has been 
used. The time evolution of Pj obeys dpj/dt — [pj, H]/ih. 
The equations of motion which are coupled to each other 
in the classical mechanics become decoupled in the quan- 
tum case where the anti-commutation relations of f play 
a key role. For example, in the conventional spin Hall 
system, the external electric field is applied to drive the 
spin current and dpi/dt = eEi. The equation of motion 
for f reduces to the one ever considered in Ref. [f| . 

In the above, we investigated the motion of a particle 
with spin in the Yang-Mills field from a classical point of 
view which differs from that of Ref ^(|, and presented 
a much more clear picture for the non-conservation of 
spin current, which has not been correctly exposed be- 
fore. We revisited to the classical counterpart of continu- 
itylike equations for spin-current density, which has the 



same form as we ever proposed quantum-mechanically 
in a previous paper §)]. With those results, we eluci- 
dated the equations of motion for a particle with spin in 
the presence of Yang-Mills fields as well as the Maxwell 
fields. Although various authors [j| have considered the 
equation of motion, the complete form J2J has not been 
recognized yet. It is worthwhile to note that the appear- 
ance of the coupling between the spin orientation and 
the Yang-Mills fields brings about the time-dependent 
effective fields even when the Yang-Mills fields are time- 
independent. In other words, the motion of the particle 
is greatly affected by the procession of its spin. From one 
of equations of motion, we can easily obtain the condition 
when the spin relaxation time is infinite and a special case 
of our conclusion recovers the previous result discussed 
by other authors. We have taken the spin Hall system as 
a concrete case and found that the tip of the spin vector 
undergoes a cyclotron motion with the frequency deter- 
mined by the strength of the spin-orbit coupling and the 
initial velocity of the particle. Furthermore, we reviewed 
the quantum counterpart of the equations of motion and 
indicated that the anti-commutation relation of the Pauli 
matrices makes the equations of motion decoupled, which 
is thus easy to solve. 

The work was supported by NSFC grant No. 10225419. 
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